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A treatment of line opacity is one of the most crucial problems for the light curve (LC) 
modeling of Type la supernovae (SNe la). Spectral lines are the main source of opacity inside 
SN la ejecta from ultraviolet through infrared range. A lot of work has been done on this 
subject [1, 3, 6, 7, 8, 9]. 

The problem itself can be divided into two parts: 

• How flux and energy equations are changed in the expanding medium; 

• How one should average flux and intensity, as well as extinction and absorption coeffi- 
cients, over energy bins used in calculations. 

Almost all previous research was about the flux equation. Here we will focus on the 
energy equation. We will suppose a free expansion of gas, i.e. linear law for the velocity 
distribution v = r/t. The principal difference between the behavior of flux and intensity in 
comparison to the static case is that the flux always becomes lower in the expanding media, 
while the sign of the change of intensity depends on the temperature gradient. It is clear 
qualitatively. When we are sitting inside the gas and would like to measure a flux and an 
intensity at a fixed frequency in the static case we should just measure how the surrounding 
gas emits and absorbs at this frequency. If this specific frequency corresponds to a strong line 
we will not be able to see distant layers of gas, since rather large optical depth accumulates 
quite close to the observer. In the continuum we can see rather distant gas layers. So at the 
frequency of a strong line we observe almost local gas without temperature gradient, hence 
we measure rather low flux value, and the intensity corresponds to the local blackbody value. 
In the continuum the flux grows appreciably because of growing gradients, and the intensity 
corresponds to the blackbody at r ~ 1, where conditions differ from local ones. 

In the expanding case, a photon's mean free path becomes less than in the static medium, 
even at the frequency which corresponds to the continuum at the rest observer frame. If we 
encounter a strong line blueward, before r ~ 1 is reached in continuum, then, most probably, 
we will not be able to see further out. So the temperature gradient at the mean free path of 
a photon becomes less than in the static case, therefore the flux drops down. The intensity 
in the expanding medium would then correspond to the blackbody radiation at the distance 
where the line forms, but it is redshifted by Doppler effect. For a constant temperature 
and pure continuum, the intensity at any space point and any frequency is lower than the 
local blackbody intensity. A strong nearby line enhances the value of the observed intensity 
(opposite to the behavior of the flux, which drops down with the comparison to a pure 
continuum case) and diminishes the difference between this intensity and a local blackbody. 
Fig. 1 shows all these dependencies for the zero moment of an occupation number, which 
is proportional to an angle averaged intensity (the zero moment of intensity): J„ = ^p-fo] 
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Figure 1: Left panel: the zero moment of an occupation number fa. Dots correspond to the 
local blackbody value, dashes show the blackbody at the distance where continuum forms, 
solid line demonstrates the behavior of /o when five spectral lines are present in the frequency 
interval we are interested in. Right panel: logarithms of monochromatic a (thick lines, "top 
hat" profile), b—j, and a(b — j). The main problem is to integrate the last function, especially 
when the number of lines is not 5, but several hundreds. 



To derive these dependencies formally, we have to solve the Boltzmann equation in the 
comoving frame for a spherically-symmetrical flow: 

idf df i-^ df udf 

"5 = Vf ~ Xff , (1) 

c ot or r an ct ov 

where we put r = vt, which is quickly established in SN ejecta. When we apply a diffusion 
limit, so that the space derivatives are negligible, the solution of the Boltzmann equation is [3] 



rtv/to v \ [ v v - 

fo= ot (tv/v, v) b(tv/u, v) exp — / a [tv / v , v) dv 



dv , (2) 



where 6 is a blackbody occupation number, and the dimensionless absorption coefficient a* = 
act. 

Last equation gives us an exact solution for a monochromatic value of /o- Numerical codes 
which are used for computations of SNe radiative evolution operate with averaged numbers, 
and each energy bin contains hundreds of spectral lines. During such computations one needs 
to solve equations for energy exchange between radiation and gas in each bin: 

d < ^ > cc < a u > (<K > - < j u >). (3) 

The problem here is to average a correctly inside a bin with hundreds (or thousands) of lines. 
The right way is to average it so that 



(4) 
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Figure 2: The distribution of J/B inside the ejecta at different times demonstrates how matter 
and radiation are decoupled. 

The main difficulty here is to derive an average value of < a(b — fo) >. As we said above, 
in the strong lines fo = b. If an optical depth in the continuum is also very high, so that 
radiation is tightly coupled with matter and fo = b, then b — fo ~ and it is not so important 
how we treat a. It becomes important when radiation is decoupled from matter. As it is seen 
in the fig. 2, they are decoupled within the entire SN la ejecta even before maximum light 
(which occurs about 20th day after explosion). A similar result was obtained in [5]. So the 
correct treatment of < a > is essential. We need to take the integral (2) accurately. It is very 
hard to do numerically, since a jumps up and down by many orders of magnitude hundreds 
of times. Moreover, where a is high, there b — fo is extremely small (see fig. 1). Analytically, 
one can solve this problem, but only for very simple line profiles. We assume that lines have 
a rectangle shape ("top hat" profile). For rectangle lines in the Rayleigh-Jeans regime and 
infinite expanding medium (so that at given v one can observe emission with arbitrary shifted 
Ai/ in the transparent case) we have 
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To get < a > from (4), we need to integrate this expression once more. This also can be 
done analytically, and adds the third summation over all lines in the bin. Computing of these 
long sums is very time consuming. Probably, the efficiency can be improved along the way 
proposed in [2]. 

M=0.5M o t=20d T=10.000K M=0.5M a t=100d T=5,000K 
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Figure 3: Logarithm of absorption coefficient calculated accordingly to our formulae (solid) 
and in Eastman-Pinto approximation (dotted). Top panels show the absorption for two sets 
of (time, temperature) inside the layer of W7 model which consists of mostly 56 Ni, while 
bottom panels correspond to the layer of intermediate mass elements. 

Using our expression for < a >, we have computed the opacity tables for standard SN la 
model W7. The results are compared with the absorption derived with the Eastman-Pinto 
approximation [7] (see fig. 3). In the UV region, where free-free and bound-free emission 
dominates, < a > remains almost unchanged, while it differs up to a factor of two in the 
range between near UV through IR. 

LCs for W7 (fig. 4) were calculated with the help of the code STELLA [4]. UBV LCs 
remain almost unchanged, while I-band is affected appreciably (fig. 4). This is explained by 
the changes of opacity and spectrum in the I-band. 
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Figure 4: UBVI light curves of the model W7 with two opacity approximations: Eastman- 
Pinto (black dashed) and the way of averaging presented here (red solid) . UBV bands remain 
almost unaffected, while I-band shows now a complicated structure, more similar to what 
one can observe for real SNe la, though not identical. Crosses, stars and triangles show the 
observational LCs for three SNe la. 
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